Implementing Lowest-Order Methods for Diffusive Problems with a DSEL

Industrial simulation software have to manage: (i) the complexity of the underlying physical models, (ii) the complexity of numerical methods used to solve the PDE systems, and finally (iii) the complexity of the low level computer science services required to have efficient software on modern hardware. Nowadays, some frameworks offer a number of advanced tools to deal with the complexity related to parallelism in a transparent way. However, high level complexity related to discretization methods and physical models lack of tools to help physicists to develop complex applications. Generative programming and domain-specific languages (DSL) are key technologies allowing to write code with a high level expressive language and take advantage of the efficiency of generated code for low level services. Their application to Scientific Computing has been up to now limited to Finite Element (FE) methods and Galerkin methods, for which a unified mathematical framework has been existing for a long time (see projects like Freefem++, Getdp, Getfem++, Sundance, Feel++ [3], Fenics project). In reservoir and basin modeling, lowest order methods are promising methods allowing to handle general meshes. Extending Finite Volume (FV) methods, Aavatsmark, Barkve, Bøe and Mannseth propose consistent schemes for non orthogonal meshes while stability problems are solved with the Mimetic Finite Difference method (MFD) and the Mixte/Hybrid Finite Volume methods (MHFV) [1]. However the lack of a unified mathematical frame was a serious limit to the extension all of these methods to a large variety of problems. In [2], the authors propose a unified way to express FV multi-points scheme and DFM/VFMH methods. This mathematical frame allows us to extend the DSL used for FE and Galerkin methods to lowest order methods. We focus then on the capability of such language to allow the description and the resolution of various and complex problems with different lowest-order methods. We validate the design of the DSL that we have embedded in C++, on the implementation of several academic problems. We present some convergence results and J.-M. Gratien IFP Energies nouvelles, 1 et 4 av Bois Préau 92500 Rueil-Malmaison, France e-mail: [email protected]